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In this article we demonstrate how narrow stripes and pairing follow from a single simple principle: 
a strong antiferromagnetic background forces injected holes to hop in steps of two such that they 
always remain on the same sublattice. When applied to a domain wall in an antiferromagnet, 
this simple effect naturally leads to a number of results. First, we demonstrate that the holes are 
immediately (exponentially) localized on stripes. Consequently, we show that the holes on a stripe 
favor the formation of pairs on neighboring rungs or sites. We recover in this way the well known 
pictures of stripes found by White and Scalapino by numerical Density Matrix Renormalization 
Group (DMRG) calculations. Throughout this work much emphasis is placed on the problem of a 
two leg ladder immersed in a staggered magnetic field. 

PACS numbers: 71. g 



I. INTRODUCTION 

It has been more than 15 years since the discovery of 
the high-T c superconductors^. Although much has been 
learned since, there is still no satisfactory explanation of 
what causes the superconductivity. There is widespread 
agreement that it should be possible to describe the elec- 
tronic properties of the square Cu02 lattices by the two- 
dimensional repulsive Hubbard model. This well known 
model is given by 



Scalapino^ have also shown that the ground-state for 
stripes in the t-J model is approximately quarter-filled. 

The question of why stripes are quarter filled is not, 
however, the only enigma surrounding stripes. Other 
questions that still need to be answered are for instancei: 
why they are so narrow? Are they insulators or metals? 
Are they bond-centered of site-centered? Are they super- 
conducting? Is the conductivity perpendicular or along a 
stripe? Are they like Akp or like 2k p density wavesZ? In 
an offshoot to this article we will detail a possible reason 
why stripes are quarter-filled. 
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where c\ a creates an electron on site i with spin a and j 
is a nearest neighbor of i. This model contains both the 
movement of the electrons (hopping) (t, kinetic energy) 
and the interactions of the electrons if they are on the 
same site (U, potential energy). The Hubbard model is 
one of the simplest possible models of interacting elec- 
trons. The main problem is then finding the possible 
solutions of this model. This however has proven not to 
be an easy task. The Hilbert space of the model is too 
large for exact numerical solutions. Monte Carlo simula- 
tions have the minus sign problem. U is approximately 
equal to 8t, so we are neither in an extremely strong, 
nor in a weakly interacting limit. Therefore, perturba- 
tion theory is of limited use. In order to get somewhere 
we therefore need to invoke simplifications. In this arti- 
cle the main assumption is that we assume that a lightly 
doped system displays sublattice order. 

Stripes were first detected in cuprates^ in the famous 
neutron scattering experiment of Tranquada and cowork- 
ers. Earlier theoretical approaches already predicted 
stripes before their experimental discovery^. We know 
experimentally that stripes are quarter-filled from a plot 
of incommensuration versus doping (the Yamada plot), 
which is a straight line from to 0.125 doping with a 
slope of l/4£. The DMRG calculations by White and 



II. OUTLINE 

Throughout this work we will focus our discussion on 
the Hubbard model. Nevertheless, all that we will state 
can be applied to a wider range of models with additional 
longer range interactions and hoppings. We will assume 
from the outset that bond-centered stripes (stripes with 
the geometry of two leg ladders) form antiphase domain 
walls in the surrounding antiferromagnetic and demon- 
strate how self consistent localized eigenstates and pair- 
ing follow from this assumption. 

In Section l|III|) . we introduce the sublattice principle 
which will form the backbone of our analysis. In order 
to streamline the quintessential physics, we will be brief 
in exposing this principle. For further details regarding 
the underpinnings of this principle, the interested reader 
may peruse Appendix A. Invoking the sublattice princi- 
ple, we will show in Section (fVf> how exponentially local- 
ized wavefunctions will be found when we assume, self 
consistently, that stripes form antiphase domain walls. 
We will find that the transverse stripe scale is of the order 
of the lattice constant. What drives stripe formation in 
our picture are not confining magnetic string potentials, 
but a rather novel kinetic effect which we term dynamical 
confinement. 

Once we establish, self consistently, that bond centered 
stripes form domain walls in the surrounding antiferro- 
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magnet, we will move in small steps towards examining 
further microscopies. As we will explain in Section (|V1|I . 
we will consider a bond centered stripe engulfed by a 
surrounding antiferromagnetic region as a two leg ladder 
immersed in an external staggered magnetic field. We 
will then consider the problem of a single electron on a 
staggered empty two leg ladder (Section l|VH|) ') which will 
facilitate the analysis of a single hole on an otherwise full 
staggered two leg ladder (Section IjVlllfl). Both problems 
will lead to similar results. 

This will be followed by a similar analysis for a two 
electron system on an empty staggered ladder and the 
inverted problem of two holes on a full staggered ladder, 
in Sections (|IX(I and (Jxj respectively. The surprising 
conclusion is that in this case there is an essential differ- 
ence between electrons and holes. Numerically, we find 
that for holes, pair states are slightly favored over single 
hole states although the correlations are very faint. 

Having shown how narrow bond centered stripes with 
bound pair states emerge, in Section l|XH|) we fuse all of 
our finding together to reconstruct the real space DMRG 
pictures of White and Scalapino. 

In Section (|X111|I , we will examine the effects of longer 
range Coulomb effects and additional longer range kinetic 
terms to show how much of our self consistent analysis 
can easily be fortified by the addition of such terms. 



III. THE SUBLATTICE PARITY PRINCIPLE 

The physics of holes and spins in two dimensions is a 
very complicated subject. In this article we will invoke 
the simplifying assumption that much of the low energy 
physics of holes in a strong antiferromagnetic background 
can be summed in a nutshell: 



Holes can only move in steps of two. 



This principle is equivalent to the omission of spin flips, 
waves, and string states. This sublattice parity principle 
is, of course, a gross oversimplification- the physics of hole 
motion in an antiferromagnetic is a fascinating and rich 
topic. Nonetheless, the low energy single hole dispersion 
curves (requiring careful extensive numerical work) coin- 
cide very well with those immediately following from this 
principle. For a review of this often overlooked principle, 
the reader is invited to read Appendix 

Stated in formal terms, the principle amounts to or- 
der, which has lately been much discussed in the context 
of its viable destruction 8 . 

Many beautiful works exist on the subject of stripe for- 
mation and dynamics^. In recent years, hole dynamics 
on the stripe has been primarily addressed in terms of 
spinon-holon excitations, as in the early works of Tch- 
ernyshyov and Pryadko^i. In this article we do not as- 
sume that the elementary excitations are spinons and 
holons. Our approach is closer in spirit to the very inter- 
esting work of Chernyshev, White, and Castro NetoiS, 



which aim to answer the same fundamental questions 
concerning stripes. Their work examines the elementary 
excitations within the framework of the t-J z model by 
examining retraceable paths of a hole out of the stripe 
into the surrounding antiferromagnet. We arrive at much 
of the same physics using the sublattice parity princi- 
ple as our only guide. Our derivations are much more 
pedestrian yet physically transparent as compared to the 
detailed Green's function and numerical DMRG analysis 
carried out by these authors. 

Many works place much emphasis on the string states 
created by hole motion or on kinks 9 . In our low energy 
analysis, there are no string states. For a discussion on 
how string states can be avoided and for a demonstra- 
tion that they do not dominate the low energy physics, 
the reader is referred to Appendix ljA"}> . Similarly, the 
role of a domain wall in an antiferromagnet as an ef- 
fective attractive potential for holes simply by the bad 
ferromagnetic magnetic bonds that they remove, and the 
careful interplay between string states and transverse and 
longitudinal hole kinetics is not what we consider here. A 
large body of literature complements our simple asymp- 
totic low energy analysis. 



IV. A BOND CENTERED STRIPE IN THE 
ANTIFERROMAGNET 

Fig. OJ shows the system which we will be looking at. 
We will examine the anatomy of a bond-centered stripe 
in an antiferromagnetic background. As fortold, we will 
neglect spin-flips and assume the background to be a per- 
fect Neel antiferromagnet. We will further assume that 
there is a phase-shift in the staggered antiferromagnetic 
order parameter as we traverse the stripe. This implies 
that there is a seam of ferromagnetic bonds between spins 
on the two different legs of the stripe. The ferromagnetic 
bonds will cost a lot of energy. However, introducing 
holes into the stripe will reduce this strain. The ferro- 
magnetic bonds will turn out to be essential to explain 
the stability of stripes against "hole evaporation" . 

In order to understand why a system would want to 
form such a stripe, we make the theoretical assumption 
that the system starts out with an antiphase boundary 
as shown in Fig. |2 Theoretically, one can force such a 
system by imposing suitable boundary conditions. If the 
system is L sites long, this will entail an energy penalty 
of, approximately, 2LJ Z . Introducing holes into this sys- 
tem will reduce the energy of this ferromagnetic seam 
and ameliorate matters. 

Topologically, stripes are not literally domain walls in 
an antiferromagnet (whose spins S are continuous vari- 
ables and not Ising like) but rather topological excitations 
of the antiferromagnet known as skyrmions. As noted by 
Wilczek and Zee^ and later incorporated by others in 
rather novel theories, e.g. Wiegmann^, skyrmions in an 
antiferromagnet are cylindrical domains separating Neel 
states shifted by half a period. We note that topolog- 
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FIG. 1: Schematic representation of a quarter-filled bond- 
centered stripe in an antiferromagnet. We assume there are 
no spin-flips in the antiferromagnet. There is a ir- phase shift 
in the staggered magnetization over the stripe. The cartoon 
above is not to be taken too literally. In reality, the pairs are 
smeared along the rungs. 
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FIG. 2: We will make the assumption that the system starts 
out as antiferromagnetic with a "domain waW - a skyrmion. 
We will analyze what happens if we add holes to this system. 

ically, stripes are identical to skyrmions stretched out 
to form domain walls spanning the entire lattice. Berry 
phase effects in the antiferromagnet associated with do- 
main wall stripes can give rise to exotic statistics similar 
to that of skyrmions in 2+1 dimensions- as suggested by 
extending the results o^ to our domain walls. 



V. DYNAMICAL CONFINEMENT (ON-STRIPE 
HYPOTHESIS) 

We now examine what happens if we introduce a 
hole into the antiferromagnet with a domain wall (see 
Fig. 0). Once again, our major assumption still is 
that this hole moves in steps of two such that it will 
always stay on the same sublattice. We reiterate a 
note for the experts: our analysis has nothing to do 
with string states and kink (wiggle) motion beauti- 
fuly analyzed by manjiiSiil, and many of the works 
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parture -the sublattice parity principle- is just a way to 
explicitly avoid high energy string states from the out- 
set). Nor, does our bare analysis have anything to do 
with the role of the magnetic alleviation energies for a 
hole drifting into the stripe (whereby it removes bad fer- 
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FIG. 3: A number of holes and an antiferromagnetic domain 
wall immersed in the surrounding antiferromagnet. We will 
show that the holes will automatically end up on the domain 
wall because of the single step up/down movement with an 
amplitude t which is possible within the domain wall (the 
stripe), but not in the surrounding antiferromagnet. 



romagnetic bonds) which plays the role of an effective 
attractive potential. Such effects will only enhance the 
bare findings that we report here. 

With these remarks in mind, we note that if the hole 
is in the antiferromagnet, far away from the stripe, it has 
eight positions were it can move to (see Fig. |3J). If it is 
sitting next to the boundary (just outside the stripe), it 
loses one hop, because of the reversed sublattice struc- 
ture. The remaining seven hops are all the same as if the 
boundary was not there. 

Larger changes happen with a hole located within the 
stripe. For such a hole there are only six hops left. In 
Fig. |21 the hole on the stripe cannot move downward left 
or downward right because of the change in sublattice 
structure. So, naively, one might come to the conclusion 
that a hole has more kinetic energy in the antiferromag- 
net then if it is on or near the boundary/stripe. Naively, 
confinement of the hole in the stripe would cost kinetic 
energy. However, this is not the case. 

The reason is that within the antiferromagnet the holes 
can only move in two steps longitudinally with an am- 
plitude to or diagonally with an amplitude £2 which are 
of the order of J = t 2 /U <C t (the uniniated reader is 
refered to Appendix^]) for an exposure to the origin of 
the exchange coupling J). On the stripe the hole can 
stay on the same sublattice by doing a single (instead of 
double) step up or down with amplitude t. As a conse- 
quence, the hole has a much larger kinetic energy if it is 
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situated on the stripe. Therefore we find that the hole 
automatically drifts toward the stripe: the amplitude for 
the ground state wavefunction is maximal for sites on the 
stripe and decays exponentially the further the sites are 
away from the strip (see Fig. 0}. Holes are driven by an 
increase in kinetic energy (not exchange energy) to the 
domain wall. We term this mechanism "dynamical con- 
finement . The word dynamical is used because it is the 
motion of the holes on the stripe that lowers the energy. 
The primary role of kinetic energy in favoring the stripe 
as a ground state has been discussed many times before, 
e.giiSiiLiS. Nonetheless, the trivial (yet sizable) lowering 
of the kinetic energy allowed by hole motion along the 
ferromagnetic rungs on the domain wall seems to have 
gone unnoticed. 

For the two dimensional problem of a single hole in the 
vicinity of a full domain wall embedded in a surround- 
ing antifcrromagnet (of the variant shown in Fig. JHl), 
we may work in a very much reduced Hilbert space. In 
effect, the strongly correlated many particle system re- 
duces to a single particle problem defined on half of the 
lattice. As there are, ab initio, four possible electronic 
states on each site (empty, one electron with its spin up, 
an electron with spin down, and an up and down elec- 
tronic pair), the Hilbert space of an N site system is, 
trivially, of dimension 4 . Invoking the sublattice parity 
principle allows us to reduce the Hilbert space to a mere 
size of N/2. The Hubbard Hamiltonian is then reduced 
to a purely kinetic model, having an amplitude £ for near- 
est neighbor direct hops within the stripe, amplitudes to 
and £2 f° r longitudinal and diagonal two step hops (as 
depicted in Fig. [3] for the horizontal motions) within the 
antifcrromagnet and an amplitude t\ for two step hop- 
ping within the stripe along the axis of the ladder. 

Before doing a more detailed analysis, let us first give 
the reader an intuitive feeling of where we are heading. 
Looking at Fig. ©, we note that a hole in the bulk 
disperses with an energy 

e bulk = eo - 2£ (cos2fc x + cos2fc y ) 

-2£ 2 ((cos(fc :E + k y ) + cos(k x - ky)). (2) 

Similarly, along the ladder, 

e stripe = eo ± * - 2£i cos 2k x , (3) 

where the second, ±£, contribution denotes the ener- 
gies of the bonding/antibonding states along the rung. 
Within this approximation, a sizable finite gap A w £ 
separates the minima of both dispersions. In our analy- 
sis, it is this 0(t) gap which is the main driving force for 
self consistent stripe formation having nothing direct to 
do with string states and magnetic alleviation energies 
(which are all of order O(J)). 

In reality, the hole can hop between the stripe and its 
surrounding bulk: the stripe states may be connected to 
the bulk. The low lying stripe states are however much 
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FIG. 4: Groundstate wavefunction perpendicular to a 
domain-wall for a single hole in an antiferromagnet. This 
was calculated using the simple approximation discussed in 
the maintext. The hole prefers to be on the stripe. We find 
H'\ ~ I ~f \ v ■ The line in the figure just connects the midpoints. 
N x and N y are the linear extents of the two dimensional sys- 
tem and the hopping amplitudes to, t\ and £2 are as defined 
in Fig. (3). 

lower in energy than their counterparts within the bulk. 
As a consequence, holes become localized on the stripe. 
As t ^> t 2 /U = O(to) = 0(t<z), the on stripe dispersion is 
markedly lower in energy than its counterpart for motion 
within the antiferromagnet. This leads to an exponential 
decay of the lowest eigenstates states out of the stripe. 

We have numerically diagonalized a 12 x 21 system 
with a direct nearest neighbor hopping amplitude 4=1, 
effective longitudinal and diagonal next nearest neighbor 
hopping amplitudes to = £2 = —0.1 within the antiferro- 
magnet, and a next nearest neighbor amplitude along the 
axis of the stripe (ladder) t\ = —0.1. The groundstate 
wavefunction is, to good numerical accuracy, given by 

V^,2/)«|y|V fc ^. (4) 

Fig. (@J depicts this wavefunction along the direction 
perpendicular to the stripe. This figure is similar to Fig. 
22 from Chernyshev et. <ilM* who calculated this from 
a 11 x 7 t-J z system employing the numerical DMRG 
method. Our approach leads directly to this exponential 
confinement due to the dynamical confinement, without 
the need for complicated numerical calculations. 

We will now aim to give the reader an intuitive grip on 
the physics along with the ability to easily derive rigorous 
bounds on hole localization within the planar problem. 
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If we were to allow a hole on the particle to hop to 
a larger number of neighbors within the bulk (i.e. to 
artificially enhance kinetic motions and lower confining 
tendencies), then we may map our system into the dy- 
namics of pseudo-spin 1/2 particles in the presence of 
a transverse magnetic field. The mapping is as follows: 
let us mark each point within the bulk by its sublattice 
parity (up/down) number as depicted in Fig.@. Let us 
now tile the plane into vertical 2x1 domino blocks lying 
with their long side parallel to the y axis. Each domino 
within the bulk contains one up and one down site; the 
labelling of the up and down sites on the dominos lying 
along the rungs of the ladders may be done arbitrarily. 
Next, let us envision replacing each domino by a single 
spin-1/2 particle: the number of the fictitious spin-1/2 
particles is equal to half of the number of sites which we 
now label by m and n. For simplicity we will now assume 
that all second order hops are of equal magnitude a. Let 
us now consider the Hamiltonian 

H=-a ^ ( c lnd C nd + ci d C md )- ^ hexfDn, (5) 
{mn) ,d n£ladder 

with the fictitious external magnetic field h ext = ta x e x 
oriented along the transverse x-axis, and D n the spin of 
the psuedo-particle (domino) along the n-th rung od the 
stripe and d is the pseudo-spin polarization label. In the 
first term m and n span the entire plane; the sum is per- 
formed over all nearest neighbor site (mn) and on four 
of the eight next nearest neighbor sites (such that all 
same sublattice hoppings of the holes within the bulk are 
accounted for). For holes not far from the stripe or on 
it, the Hamiltonian of Eqn.JSJ introduces additional un- 
physical motions. These allow hoppings of the holes off 
and on the stripes which are disallowed- such additional 
terms can only enhance derealization tendencies. The 
second term in Eqn.(|3J codes for the nearest neighbor 
up/down hops along the same rung- these alter a pseudo- 
spin up state to a down state and vice versa. We trivially 
observe that the low energy dynamics corresponds to the 
motion of a | — ») particle polarized along the applied 
transverse field sensing an attractive confining potential 
of strength t along the stripe. In physical terms, the 
psuedospin polarized state | — ►) = 2~ 1 / 2 (| f) + | J.)) cor- 
responds to the symmetric low energy bonding state (and 
its counterpart | <— } corresponds to the high energy an- 
tibonding state). The ground state is simply that of a 
polarized | — >} particle (the symmetric bonding state) 
sensing a well of infinite extent along one axis and hav- 
ing a finite extent (2a, with a the lattice unit) along the 
transverse direction. As the problem is translationally 
invariant along the ladder axis, the longitudinal momen- 
tum k x is a good quantum number. The lowest lying 
wavefunction is translational invariant along the infinite 
cylindrical axis (k x = 0) and, at long distrances, has 
the transverse profile of a particle of an effective mass 
m e ff = 1/(2ct) subject to the influence of a potential well 
of depth t and width (2a) . This trivially leads to an expo- 
nentially decaying amplitude in the direction transvserse 



to the stripe. As t 3> J, this confining tendency is much 
more profound (and physical) than most of the common 
magnetic bond (J) arguments prevalent in the litera- 
ture. The reader should bare in mind that our point of 
departure- the sublattice parity principle is correct only 
as low energy scales (as compared to J). Nonetheless, 
what drives localization in our picture at low energy scale 
are the much more significant kinetic effects which out- 
shadow the common bad bond (J) counting arguments. 

The solution to this potential well problem is stan- 
dard. Scaling back by a factor of two to the original 
(non-domino) coordinates along the y~ axis, we find, 
asymptotically, 

ip(x,y) wAexp[-2a|y|]; \y\ > a, (6) 

with the lowest lying bound state of Eqn.JSJ is obtained 
by the solution of the transcendental equations 

[(3a tan (3a = aa] and [a 2 (a 2 + (3 2 ) = —t] 

a 

or [(3a cot (3a = -aa] and [a 2 (a 2 + (3 2 ) = -t] (7) 

a 

having the largest value of a. 

Eqn.© gives an upper bound on the leakage of the 
"hole" wavefunction out of the stripe. If we set a = max 
{to,ti,t2}, we obtain a rigorous upper bound on \ip\, at 
large distances \y\, consistent with our numerical findings 
(Eqn.Q). As the Hamiltonian of Eqn.© contains ad- 
ditional unphysical hopping processes from the ladder to 
its environment, a localized state found for the Hamilto- 
nian of Eqn.JSJ) implies even more localized states for the 
more restricted physical problem. 

We may similarly address the problem of an arbitrary 
stripe configuration. In general, adding additional hop- 
pings transforms the problem of the hole motion within 
and/or near stripes to a kinetic planar problem of a hole 
coupled by Zeeman couplings to magnetic fields piercing 
the plane along the stripe trajectories. A multitude of 
viable self consistent minima of narrow stripes of various 
geometries are found. Stripe dynamics is akin to the mo- 
tion of the fields (solenoids) piercing the plane; these, in 
turn affect hole dynamics by a Zeeman like effect. The 
evolution of a hole-stripe system may be addressed via 
such a self consistent scheme. 

A related way of immediately deriving lower bounds 
on the diffusion of the hole out of the stripe amounts to 
looking at the transverse cross-section of the stripe as 
depicted in Fig. (JjjJ, and examining single hole motions. 
Solving the Schrodinger equation for this one dimensional 
system, we immediately obtain a localized bonding state 
of the lowest energy. As here fewer hops are accounted 
for than in the original physical problem, the solution 
serves as a lower bound on "hole" wavefunction leakage 
out of the stripe. When fused with the previous upper 
bounds from the pseudo-spin system, we immediately ob- 
tain both lower and upper exponential bounds on asymp- 
totic hole propagation out of the stripe. 
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FIG. 5: One-dimensional example showing the reason for dy- 
namical confinement. Outside the stripe, the hole moves in 
steps of two with an amplitude to. Inside the stripe it can 
make a direct hop with amplitude t with t 3> to- 



In this scheme, splintering the sites according to their 
sublattice numbers (equivalent to the magnetic field 
sensed by a spin of definite polarization if it is placed 
at various locations throughout the lattice) introduces 
hopping parameters for the low energy dynamics with no 
potential energy in sight. As shown by our analysis, mag- 
netic alleviation effects are not imperative for achieving 
this dynamical confinement (even when direct magnetic 
effects are removed we nunmerically attain exponential 
localization with little noticable change). We retiterate 
that we are examining the reduced Hilbert space where 
no magnetic string states exist from the outset. 

In our approach, we first create an antiferromagnctic 
spin structure with domain walls and then dope the sys- 
tem. Then, if the spin structure is strong enough, the 
holes will migrate to the walls and will form a charge 
structure. Hopefully we have convinced the reader that 
sublattice parity order can indeed drive stripe formation. 
The energy scale associated with this discrete sublattice 
parity (Z2) order— can indeed be quite high: the persis- 
tence of incommensurate peaks up to at least the high 
energy resonance peak in YBCO might be interpreted as 
an indication that the stripe presits to be a domain wall 
up to very high energy scales 1 ^. Experimentally, charge 
order is a far greater robust driving force for stripe for- 
mation than spin order. The role of sublattice parity 
(albeit its high energy scales) is not clear at the time of 
writing. The reader should consider our argument as a 
self-consistent one in which the creation of the no-hole 
domain-wall first is a theoretical device that simplifies 
theoretical discussion. Summarizing, we have seen that 
once there is a ferromagnetic seam in a very strong anti- 
ferromagnet, single holes will automatically move to this 
seam and be exponentially localized onto them. In the 
remaining part of this article we will see what the conse- 
quences are of this result for more than one hole on the 
stripe. 



VI. STAGGERED LADDER SYSTEMS 

We have seen that by dynamical confinement we may 
consider an antiferromagnetic domain wall as a two-leg 
ladder with staggered boundary conditions. So, once 
again we can restrict our Hilbert space. Ladder sys- 
tems have been extensively studied throughout the years. 
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FIG. 6: Single electron on an empty ladder. In principle the 
electron can hop to all its nearest neighbors with an amplitude 
t. However, the ladder is embedded in an antiferromagnetic 
background that influences the electron on the ladder. 



Hundreds of works on standard (unstaggered) ladder sys- 
tems have been carried. For a well-known review see 
Dagotto and Rice 40 . 

In our case, the influence of the antiferromagnet sur- 
rounding the ladder is still there. The surrounding anti- 
ferromagnet effectively gives rise to staggered boundary 
conditions. In the up and coming, we will look at two- 
leg ladders with staggered fields mimicking the surround- 
ing antiferromagnet. The staggered fields endow the lad- 
der with a sublattice parity structure. Unfortunately, to 
date, staggered ladders have not been investigated inten- 
sively. 

Krotov, Lee, and Balatsky 14 examined staggered lad- 
ders for the Hubbard model at small U. In this limit, one 
starts with a non-interacting system and then derives the 
renormalization group equations. Staggered spin ladders 
without any holes have been studied by Wang et. a@&. 
However, staggered ladder systems with holes in the large 
U limit have not yet been addressed before as far as we 
know. We will now analyze staggered two leg ladders in a 
rather pedestrian manner. We verified the perturbative 
results that we will cite by explicit numerical computa- 
tions. 



VII. THE SINGLE ELECTRON ON AN EMPTY 
STAGGERED LADDER 

To understand the effect of the boundary conditions 
on the movement of electrons and holes on the ladder we 
first examine a single electron on an empty ladder im- 
mersed in a staggered magnetic field. This is, of course, 
a hypothetical situation, but it is nevertheless a good 
starting point for our discussion. If the electron is on the 
correct sublattice (i.e. if its spin polarization is opposite 
to that of the applied staggered magnetic field) then its 
magnetic energy will be eo = — J. Being on the wrong 
sublattice leads to a magnetic energy penalty of e\ = + J. 

This gives rise to staggered potentials on the ladder. 
The resulting Hamiltonian for the system as shown in 
Fig. JJJJ, reads 



7 



j+t 



H = 



f 1 
c 1 


—t 


—t 














-t 


ei 





-t 












-t 





eo 


-t 


-t 













+ 

— c 


4- 

— 1 


(() 


u 


+ 












-t 





ei 


-t 


-f 












-t 


-t 


Cl 





-t 












-t 





eo 


-t 


-i 













-f 


-t 
-t 


eo 


-i 




ei 

-/; 



-i 
-i 



(8) 



Here we express the Hamiltonian in terms of ordered 
real-space basis states where we place the upper and 
lower sites of each rung adjacent to each other. An elec- 
tron at any site can always hop to three other sites. Be- 
cause this is a single particle problem with a translational 
invariant potential (with unit cell of size two), the eigen- 
functions of this Hamiltonian can exactly be determined. 
They are exactly given by 



ip= (B,±B, Ae ik ,±Ae ik , Be 2ik ,±Be 2ik , Ae mk , ...)• (9) 
Here A and B are determined by 

Ek = e ±t — 2t—cosk 
A 



E k = e 1 ±t-2t—cosk. (10) 
B 

The energy can immediately be obtained from 




t+j 




FIG. 7: Energy diagram for a single electron on an empty 
staggered ladder. If the staggered potential is large (large J), 
the splitting J is more important than t. In the lowest state, 
there are only N/4 states. These are the same states as for 
a single one dimensional model where an electron moves in 
steps of two. If J -C t (right) we still get the same lowest 
sector. 



—J. The other half have the hole on the wrong sublat- 
tice, with an energy of + J. Because the upper and lower 
leg of the ladder are exactly equivalent, the electron will 
slosh back and forth between the upper and lower legs. 
This, in turn, leads to bonding/antibonding linear com- 
binations of the upper and lower sites along each rung, 
further splintering the Hilbert space into two additional 
subsectors. The lowest sector, which contains N/4 ba- 
sis states, consists of wavefunctions with the electron on 
the right sublattice (e.g. a spin up polarized electron on 
the up sublattice) with the symmetric linear combination 
(bonding) of upper and lower sites. 
In the large J limit, Eqn. (I12|) simplifies to 



e a ±t — E k 
—2t cosfc 



— 2t cos k 
ei±t-E k 



Ei. 



= 0, 



(11) 



-J-t-2- 



leading to 



E k = ±t±y/j 2 + (2iC0sfc) 2 

= ±t± \/j 2 + 2t 2 + 2t 2 cos2fc. 



(12) 



This dispersion illustrates that, due to the staggered 
boundary conditions, the electron effectively moves in 
steps of two. Fourier transforming the energy shows that 
because of the symmetry between k and 7r — k, only terms 
which have an even number of steps from the starting 
point may have an amplitude unequal to zero: odd, sub- 
lattice parity interchanging, hops are banned. In more 
physical terms, this trivially corresponds to the halving 
of the period (in k) in the dispersion E k . 

Eqn. I|12l) shows that, for unphysically large J ^> t 
(a strong influence of the surrounding two-dimensional 
antiferromagnet), the Hilbert space splits up in four dif- 
ferent sectors. First, there is a splitting because of the 
boundary conditions: half of the sites have the electron 
on the right sublattice, which leads to a low energy of 



2J 



2 — cos 2k. 
2J 



(13) 



Not surprisingly, we see that we this energy corre- 
sponds to the electron populating the correct sublattice 
(—J) and being smeared along the two rungs in a sym- 
metric bonding fashion (— t). We can easily read off that 
the effective hopping is in steps of two from the cos 2k 
term. Because we are in the large J limit, the hopping 
amplitude in second order perturbation theory is given 
by |j. We have to hop twice (t 2 ) over an intermediate 

.2 

state with energy 2J. The extra contribution — 2jjto 
the energy comes from virtual excitations where the spin 
moves one position to the left or right and immediately 
returns back. There are two possible ways of doing this, 
and the amplitude again is |j. Thus, we can intuitively 
understand every term in this limit. 

There is a gap of 2t separating the lowest Hubbard 
band from the antibonding states. Because 2t is a very 
large energy scale (approximately 0.5 eV, or 5000 K), we 
can neglect the influence of the other sectors and may 
restrict our attention only to the lowest sector. 

Our discussion above hinged on the assumption that 
J is very big (J ^> t). This assumption is not satisfied 
in the physically relevant region wherein J w 0.25t. In 
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FIG. 8: Effectively, the electron moves on a one-dimensional 
chain in steps of two. There are only two effective parameters: 
the on-site energy eo and the two site hopping amplitude a. 
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the physically relevant regime, we still have a splitting of 
the Hilbert space into four sectors (see Fig. @). The 
lowest Hilbert sector is still the same one as the one for 
the large J limit. For small values of J compared to t, 
the groundstate properties are still dominated by the — J, 
—t sector. However, the low lying excitations out of the 
lowest sector are now to the sector where the electron 
goes to the wrong sublattice (a single hop) instead of 
going from a bonding to an antibonding state. Such hops 
will lead to string states. Nevertheless, as long as the 
width of the lowest Hubbard band is small compared to 
2J (which itself is also a large energy, of the order of 
2400K), we can still think in terms of bonding states 
on the correct sublattice. Within the lowest Hubbard 
sector, the electron will only be on the even lattice sites. 
Effectively, only second order hops are present. Because 
of the hybridization of up and down sites on a rung, it is 
sufficient to consider the one dimensional model shown 
in Fig. ©. 

Just as for the single hole in the two dimensional anti- 
ferromagnet, discussed in Appendix (|A"|l . a lone electron 
on the empty ladder moves as an effective quasiparticle, 
with a sublattice hopping amplitude a. 



VIII. ONE HOLE ON A FULL STAGGERED 
LADDER 

We now examine one hole in an otherwise full ladder 
immersed in a staggered external magnetic field. The 
single spin problem was very simple to solve analytically 
given our assumptions. The single hole problem is ex- 
ceedingly more difficult because it is a strongly interact- 
ing many body problem. For a system of only 6x2 sites, 
there are already approximately 800, 000 states! Systems 
of this size and larger may be addressed by employing the 
Lanczos method of diagonalization, which only finds the 
lowest eigenstates of a matrix. We employed this method 
for systems of up to size 8x2. The complicated results 
that follow from this numerical study can be understood 
quite easily for large values of J. Though not the phys- 
ically relevant regime, just as in Section (VII), both for 
J > t and J < t, the ground state is in the same sector 
of Hilbert space. 

For large values of J, the basic properties are once 
again those of a single quasiparticle. From the outside, 
we again assume that the hole cannot leave the stripe. 
An electron next to the hole can move to the position of 
the hole, leading to a string-state with one wrong spin. 



FIG. 9: We assume a single hole can move on the stripe and 
cannot leave the stripe. The surrounding antiferromagnet 
leads to staggered boundary conditions. 



The hole can also move up and down. And it can move 
in steps of two to the left and right. In addition, we can 
have spin flips on the ladder. However, just as in the 
two-dimensional antiferromagnet we neglect the spin-flip 
processes. For large values of J we find once again a 
number of Hubbard sectors. The lowest sector has N/2 
states where the hole is on the right sublattice and an 
energy of approximately —2(N — 1) J. Above this low- 
est Hubbard sector we have the Hubbard sector with a 
string state of length one. There are 2N states in this 
sector. The gap between this sector and the lowest sec- 
tor above it is A = 4J w 4800K. We could ingnore 
this sector if this gap is large with respect to the internal 
splitting of the lowest Hubbard sector. This bandwidth 
is given by At" with t" — u t ^ 4J -C t. For U ~ 8t we have 
t" k O.lt: both J and t have approximately the same 
amplitude. However, for the groundstate properties, the 
relevant properties are those of the lowest a single hole 
moving in steps of two. The band, shifted down by t 2 /2J 
by virtual hopping, has the dispersion of an inverted co- 
sine. The gap A = t, separating the on-stripe states from 
those in the surrounding antiferromagnet, is large with 
respect to t" . The spins adjust to surrounding antiferro- 
magnetic. In the [/, J> 1 limit, 



eo = -J-t-2-, 



(7 = 2, 



t 2 



U + 2J 



(14) 



The motion of a single hole on a stripe and its motion 
in a two dimensional antiferromagnet is identical. The 
only difference between the parameter sets, is sparked by 
the presence of the direct hopping term t allowed within 
the topological domain wall. This is in contrast to one 
dimensional effective theories that assume that the hole 
effectively becomes a holon on a stripe (see for instance 
Tchernyshyov and Pryadkc 1 "). Our hole is not a holon 
in the sense that it will remember if it was injected for 
a spin up or a spin down electron: they move on diffent 
sublattices. 

We see that we can consider a single hole on a stripe 
as a quasiparticle, just as a single electron on an empty 
lattice. They can both be characterized by (different) 
values for eo and a. The main difference between a hole 
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FIG. 10: For a one dimensional ladder with strong staggered 
potentials on the boundaries, a single hole effectively moves 
in steps of two. 
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on a full ladder and an electron on an empty ladder is 
that the hole moves more slowly, because it can only move 
through double occupied intermediate states, which costs 
a lot of energy. Therefore, a is small for a hole, while it 
is large for an electron on an empty lattice. Our main 
conclusion is that a single hole and a single electron are 
effectively identical. 




£o 



IX. TWO ELECTRONS ON AN EMPTY 
STAGGERED LADDER 

The equivalence between holes and electrons is no 
longer true for more than one particle. We will show that 
the low energy properties for two electrons are different 
from those of two holes. We first consider the problem 
of two electrons with opposite spin on an empty ladder. 
Various possible geometries are sketched in Fig. (|11H . 
Because only the relative distance is important, we keep 
the position of the up spin fixed. 

The two electrons can be on the same or opposite legs. 
Six amplitudes: eo, ei , £2, c, a' and r, defined in Fig. (|llfl . 
are of relevance. We notice that the hopping amplitude 
for the electrons when they are far apart (a) is much 
larger than the hopping amplitude when they are next to 
each other (r). If they are far apart, they do not notice 
each other. Two electrons have less kinetic freedom if 
they are next to each other. In order to move passed 
each other, they have to go through a doubly occupied 
state. 

In the U 3> J limit, lowest order perturbation theory 
yields 



ei = £o 
a 



-2J-4- 



2J' 



e_ 

2J : 



f 2 



U + 2J 



(15) 



(16) 



FIG. 11: Different possible configurations of two spins on a 
staggered two leg ladder. 
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The lowest Hubbard sector contains both single elec- 



FIG. 12: Hubbard sectors for two spins on an empty ladder. 
The lowest Hubbard sector contains both single electrons and 
pairs of electrons. However, the groundstate consists of single 
electrons. Pairs are only formed at a higher energy in the 
lowest Hubbard sector. 



trons and pairs of electrons. However, we have seen that 
single electrons have a larger kinetic energy than pairs. 

Therefore, the ground state consists of electrons that 
are as far apart as possible. Fig. I|12|l shows a schematic 
picture for the density of states. Only at the higher en- 
ergies of the lowest Hubbard band do we find pairs. The 
pairs have a small bandwidth because of the small hop- 
ping amplitude. Once again we find various sectors as 
shown in Fig. l|12fl . 



X. TWO HOLES IN A FULL STAGGERED 
LADDER 

To see the difference between electrons and holes, we 
now look at two holes in a filled ladder. When far apart, 
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FIG. 14: Two holes can move together without creating an 
intermediate double occupied state. Therefore, holes want to 
form real-space pairs. A hole pair can move if a single electron 
makes two hops. 
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FIG. 15: We can neglect the spins on the ladder because they 
cannot change. This is caused by the boundary conditions. 

only of the order of 4 J. Therefore, a' is much large than 
a. This leads to the formation of pairs on the stripe. 

For large U and J, employing the same convention as 
before, we find in perturbation theory, 



FIG. 13: Different configurations of two holes on a stripe 
with the relevant on site energies and hopping amplitudes. 



two holes (quasiparticles) do not notice each other. Just 
as for two electrons on an empty lattice, there are two 
states in which they strongly influence each other. The 
main question we wish to address is whether two holes 
come together and form a real-space pair, or if they tend 
to be as far apart as possible. 

If the two holes are far apart, they will move inde- 
pendently with an amplitude a = v \. 2 j ( see Fig. 1 (151 ) ■ 
However, this is different if they are sitting next to each 
other. In this case no double occupied state is needed. 

,2 

Therefore the hopping amplitude r = thus r 3> cr. 

If the two holes are on the same leg, sitting next to 
each other (see Fig. JO}), than the two holes can move 
together with a much larger amplitude than if they move 
alone. 

In the left subfigure, a spin up electron moves left. 
The intermediate state in the middle figure has a higher- 
energy of the order of J. From this position, the electron 
with spin up can move back to the starting position, or 
it can move on to end up as shown in the right figure. 
Once it is there, the hole-pair has effectively moved one 
site to the right. This figure only illustrates one possible 
hopping sequence, there are also other possible ways of 
hopping. For discussions of pairing one should be careful 
with fcrmionic minus signs sparked by the interchange of 
identical spins once a pair is made to go arouncU^iii For 
the movement of a single hole, there is an intermediate 
state with an energy of U. However, if the two holes are 
next to each other then the intermediate state will not 
have a doubly occupied site, so therefore this energy is 
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(17) 



U + 2J' 

with TV = 2L the total number of sites on the two leg 
ladder. 

The motion of the pairs on a ladder and the motion 
of single holes is coupled: both are present in the lowest 
Hubbard sector. If they form disjoint Hilbert spaces, we 
will find 75 = 62 — 2rcosk for the real space pair. For 
two single holes, E = eo — 2ercosfc, within perturbation 
theory. Because a is very small, the spectrum for single 
hole states does not exhibit much dispersion. Taking the 
energy values from Eq. I|18[). the dispersion for two holes 
on a stripe with staggered potentials is suspected to look 
like that shown in Fig. H16fl . 

The big difference between holes and electrons is that 
for two electrons the separated electrons have a large 
bandwidth and the pairs of electrons are a high energy 
state with a small dispersion. For holes this is reversed: 
the pairs have the largest bandwidth, while the separated 
electrons form a low energy excitation (string states) with 
a small dispersion. 

The lowest state for the pair cosine band is located at 
£2 — 2r and the lowest state for the separated holes is 
located at eo — 2a. In perturbation theory, both energies 
amount to 



(TV 



(2TV-8)- 



(18) 



U + 2J 

However, here we assumed the pairs and separated 
holes to be independent. But where they overlap in en- 
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FIG. 16: Theoretical dispersion relation for two holes in a 
bond-centered stripe. The groundstate consists of pairs. 



ergy, which in this case is the bottom of the band, the 
wavefunctions will hybridize. 

If we think of a hybridization of up and down sites 
along the rungs, then no possible interchange can be 
performed between the holes leading to fermionic mi- 
nus signs that elevate, rather than depress, the pairing 
energj»*. 



XI. LANCZOS CALCULATIONS 

Lowest order perturbation theory for two holes on a 
staggered ladder leads to the dispersion shown in Fig. 
(15). The ground state properties cannot be easily de- 
termined because of the overlap between separate hole 
and pair states. To test this situation numerically, we 
have performed standard Lanczos calculations for the 
Hubbard model endowed with staggered boundary condi- 
tions on a twelve site single chain with staggered bound- 
ary conditions. Solving this problem for large two leg 
ladders is inhibited by the large Hilbert space. The max- 
imum size which we are able to examine numerically is a 
6x2 system for which there are w 650, 000 states. Such a 
system is too small to observe single hole states and pairs 
accurately. Because of the strong bonding combination 
between the upper and lower leg of the ladder, the reduc- 
tion to a one dimensional line is physically justified. The 
results of the Lanczos calculation show that the ground 
state has the largest amplitude for pair states with only a 
small admixture of single hole states. Immediately above 
the ground state, the excited states are predominantly 
single hole states. There is a finite energy gap separating 




FIG. 17: If we have a quarter-filled stripe, we get the above 
picture for a stripe. 

the ground state and these excited states. Because of the 
small lattice size, we cannot determine the value of the 
energy gap AE. 

Another approach is to assume that en, £2, &, and r are 
adjustable parameters. In that case, we can perform a 
simplified calculation on a 20 x 2 ladder having twenty 
spin up, and twenty spin down sites. In the aftermath, 
this leads to a 400 x 400 matrix that can be easily diag- 
onalized. Using this model it is relatively easy to choose 
the effective parameters such that the cosine band for 
the pairs has a lower energy than the states where we 
have single holes far apart. Our numerical results for 
the Lanczos method suggest that the gap AE separating 
pair states and single particle states is sufficiently large 
for the groundstate to consist mostly of pair states. 

Nontheless, as the gap separating the single hole and 
pair states is small, the pairing physics is, a priori, sus- 
ceptible to the addition of longer range interactions and 
hopping interactions. Such effects could easily tip the 
balance between the single hole and pair states (or en- 
hance pairing tendancies). 

XII. DMRG - PUTTING ALL OF THE PIECES 
TOGETHER 

We found that if we assume that somehow the anti- 
ferromagnet has a ferromagnetic domain wall, holes will 
automatically localize on this topological line. Further- 
more we looked at the movement of a single hole on this 
staggered ladder and found that it behaves as a single 
quasiparticle that moves along the ladder in steps of two. 
Two electrons on an empty lattice try to be as far apart 
as possible. Two holes on an otherwise halfhlled ladder, 
however, form a bound pair. Fusing all of these findings 
together, we obtain the stripes as seen by the DMRG cal- 
culations of White and Scalapino^. The arguments that 
we employed throughout were very pedestrian and gen- 
eral. Throughout we invoked the sublattice parity prin- 
ciple and the related staggering potentials. Our findings 
are further supported by numerical calculations on the 
t-J model. We give a very simple physical interpretation 
for the stripes found by DMRG. 

A hole is on the stripe will delocalize by the hybridiza- 
tion of the upper and lower leg. If we place a hole on the 
upper leg, the hole will almost immediately jump to the 
lower leg and vice verse. In numerical calculations, we 
will therefore find the linear superposition of a hole and 
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FIG. 18: The hole has equal amplitude to be on the upper or 
lower leg of the ladder. Therefore, the groundstate is a linear 
superposition of a hole and a spin down. 




FIG. 19: The DMRG calculation of the t-J model by White 
and Scalapino-. This and similar calculations typically em- 
ploy open or cylindrical boundary conditions. 



the electron depicted in Fig. JTSJ) . 

White and Scalapino employed the DMRG technique 
on the t-J model 4 . Fig. 1)19(1 shows one of their well 
known results. We see that the bond centered stripes 
are very narrow. They look like a ferromagnetic seam in 
an antiferromagnet. The same result can be obtained in 
mean- field theorjii^ as shown in Fig. 1)20(1. 

The DMRG calculations of White and Scalapino 
vividly show how spin and charge nestle. White and 
Scalapino place much emphasis on the fact that again the 
holes are on next-nearest neighbor sites. On the oppo- 
site diagonal, the spins (which are on the same sublattice) 
have a very strong antiferromagnetic bonding. In effect, 
the bound pair is creating a flip in the antiferromagnetic 
lattice. This might also explain why in the exact calcu- 
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FIG. 20: Mean-field solution of the repulsive Hubbard model 
with U=8t. This calculation uses open boundary conditions. 



lations and the Monte-Carlo results the binding energy 
would decrease as a function of the lattice size. Basi- 
cally, the hole-pair is creating an antiphase boundary, 
whose energy is increasing as a function of lattice size. 
That this is indeed a good representation of stripes fol- 
lows from DMRG calculations by White and Scalapino. 
The reader is urged to focus on the central part of their 
picture- the spin and charge texture towards the bound- 
aries are more contorted by the open boundary conditions 
employed. We will assume that Fig. (|2*U|l is a useful real- 
space representation of quarter-filled stripes. We claim 
that the reason that the mean-field 1 ^ and the DMRG cal- 
culation are so similar has to do with the fact that the 
antiferromagnet is very strongly ordered. The stripes are 
so narrow because of dynamical confinement. 



XIII. LONGER RANGE KINETIC AND 
COULOMB TERMS 

Let us recapitulate the assumptions that we have in- 
voked, so far, in our analysis: 

(i) We assumed, self consistently, the existence of an 
antiphase domain wall having the geometry of a two leg 
ladder. Sublattice parity (Z?) order was assumed to pre- 
vail throughout the entire system. 

(ii) Albeit the relatively minute energy difference by 
which the lowest lying pair states were found to be fa- 
vored over single hole states, we assumed that the stripe 
was composed entirely of pairs. 

It should be noted (especially in the context of assump- 
tion (ii)), that the small energy differences we found sepa- 
rating various contending states (as well as those separat- 
ing, say, bond centered stripes from site centered stripes) 
are very susceptible to additional terms in the Hamilto- 
nian. 

The Hubbard model or its t-J truncuation is only a 
model. There are myriad very important effects that 
it does not include which could easily shift the balance 
between various nearly degenerate contending states. 

Coulomb effects (which are much greater importance 
here than elsewhere given the poor screening in these 
materials) enhance and stabilize stripe order: a uniform 
charge density order is strictly forbidden by the diver- 
gent Coulomb penalty that it will incur. This point has 
been emphasized by Emery, Kivelson, and coworkers&^S. 
Fourier transforming the Hamiltonian and looking for 
the minimizing waveumbers, we are able to see how the 
charge density modulations will evolve and lead to stripe 
order— once lattice effects are taken into account. 

Even if the stripe correlations found by the DMRG cal- 
culations of White and Scalapino on the pure t-J model 
are found to be due to open boundary condition effects, 
when long range Coulomb interactions are introduced, 
charge stripe order will be further stabilized. Given the 
natural coupling between spin and charge^, this will fur- 
ther enhance the sublattice parity flips across the stripe 
that we assumed from the outset (assumption (i)). 
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Less emphasized are the role of higher order kinetic 
terms. These can easily tip the balance: a next nearest 
neighbor hopping increases pairing significantly. Within 
the pure t-J model (with a vanishing direct diagonal hop- 
ping amplitude t' — 0) pairing correlations are infinitesi- 
mal and are further frustrated by the n phase shift across 
the domain wall. Numerically, t' which allows holes to 
move on the same sublattice enhances pairing dramati- 
cally. Furthermore, numerically, a negative t' is seen to 
favor stripe formation^. 

XIV. CONCLUSION 

We examined the consequences of hole motion in a 
strongly ordered antiferromagnet. We show that holes 
move to antiferromagnetic domain walls and effectively 
form two-leg ladders. This dynamical confinement of 
holes onto the stripe is caused by the sublattice struc- 
ture of the antiferromagnet and the increase in kinetic 
energy on a domain wall. The effective two-leg ladders 
still feel the influence of the surrounding antiferromagnet 
in terms of a staggered boundary potential. As shown by 
Krotov, Lee, and Balatsky«i4, this increases the tendency 
to superconductivity. The reason for this is that holes 
on the stripe form real-space pairs. Our results natu- 
rally lead to a theoretical description of the numerical 
results of White and Scalapino. Therefore, our start- 
ing assumption- the sublattice parity principle- and the 
neglect of string states that it entails might have some 
validity. Our article furnishes yet another example of a 
kinetic theory^- string states and magnetic (potential) 
effects of the like are not the immediate driving force for 
stripe formation and hole dynamics within our scheme. 

In a upcoming paper we will show that these real-space 
pairs can be mapped to an effective one-dimensional XXZ 
model in a transverse field. This will allow us to discuss 
the filling fractions of stripes. 

Acknowledgments 

The authors would like to thank Wim van Saarloos 
for encouragement and discussions and Jan Zaanen for 
critical reading of the manuscript and for subsequently 
pointing out to usiS anctii. We were supported by the 
Foundation of Fundamental Research on Matter (FOM), 
which is sponsored by the Netherlands Organization of 
Pure research (NWO). 

APPENDIX A: THE SUBLATTICE PARITY 
PRINCIPLE 

The sublattice parity principle amounts to the assump- 
tion that a hole in an antiferromagnet is free and can 
move to all of its next nearest but not to its direct neigh- 
bors. This principle is based on and follows from the 



sublattice structure of antiferromagnetic order. We will 
provide phenomenological proof for this assumption by 
examining the dispersion relations obtained by very de- 
tailed numerical works. 



1. Spin flips 

It is well known that the Hubbard model for a half- 
filled system in the large U limit leads to an effective 
Heisenberg model, 

^Heisenberg = ^2 ~^^ S t S j + S i s t) + J * S i S h ( A1 ) 

m 

with J_l = J z = J. 

The last (Ising) term of the Heisenberg Hamiltonian 
wants to make the spins on neighboring lattice sites point 
in opposite directions. If only this term is present, the 
groundstate will be a perfect Ising antiferromagnet with 
long-range (Neel) order: the lattice can be subdivided in 
an up sublattice and a down sublattice. However, the 
first (XY) term in the Hamiltonian can undo this order 
by flipping two neighboring, opposite spins. In principle 
this term can completely destroy the long range order 
and indeed it does so in one dimension. In dimensions 
d > 3 clear sublattice order prevails. A large amount of 
effort in numerical calculations and the very important 
analysis of the non-linear sigma model by Chakravarty, 
Halperin and Nelsonii have shown that the ground-state 
of the two-dimensional Heisenberg model (and thus of 
the undoped Hubbard model) is a long range antiferro- 
magnet for which there exists a two-sublattice structure. 
The groundstate of the antiferromagnet is not exactly 
given by the classical Neel state which one would expect 
from the Ising model. The spin-flip term leads to a fi- 
nite density of flipped spins (approximately 20%) which 
give rise to a lowering of the magnetization. It also leads 
to spin- waves. These spin- waves destroy the long-range 
antiferromagnetic order at any finite temperature. How- 
ever, one can define a correlation length and on distances 
smaller than this length, we can still speak about local 
antiferromagnetic order. For T > 0, we have local order 
with a correlation length that decreases with increasing 
temperature. Therefore, thinking in terms of an Ising 
model is not completely incorrect as long as one keeps in 
mind that there is a finite density of flipped spins in the 
lattice. In the remainder of this article we will neglect 
these spin flips. 

A moment's reflection reveals that [H, S** *] = for 
the Hubbard model (a property inherited to its descen- 
dents), and net magnetization is strictly preserved. Con- 
sequently, if Neel order prevails also when a single hole 
is introduced, hole motion is rigorosly restricted to one 
sublattice: if a hole could indeed hop to its neighboring 
site on the opposite sublattice, the magnetic quantum 
number (S z /h) would remain unaltered. On the other 
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FIG. 21: If a single hole in an antiferromagnet moves by near- 
est neighbor hopping of surrounding spins, it cannot move 
without creating a string of flipped spins. For small values of 
J_l, there are no spin flips that can destroy this string. This 
leads to a linearly increasing confining potential and the eigen- 
functions are Airy wave functions. Because of this string, a 
single hole cannot easily propagate through an antiferromag- 
net by nearest neighbor hops. 



FIG. 22: A single hole can still move through an antiferro- 
magnet without generating string states by hopping in steps 
of two. This way it stays on the same sublattice. For this it 
has to go through an intermediate state with a double occu- 
pied site. This intermediate state costs an energy U. This 
leads to movement of the hole on the same sublattice with an 
amplitude oc t 2 /U. 



2. Numerical dispersions 



hand the two low energy "vacuum" states corresponding 
to the injection of a hole on the two different sublattices 
have magnetic quantum numbers (S z /h) differing by ±1. 
These states obviously cannot be connected by any of the 
magnetization conserving processes of H. Myriad ana- 
lytical treatments of the t-J and Hubbard models were 
aware of this selection rule and have computed transi- 
tion matrix elements for a single hole between sites of the 
same sublattice. For one example amongst many see—. 

A well known argument leads to the conclusion that a 
single hole cannot propagate freely in an antiferromagnet. 
As shown in Fig. (|21|) . whenever a non backtracking 
single hole moves from one sublattice to the other, it 
creates a string of flipped spins in its wake. This leads 
to a linear confining potential that strongly inhibits hole 
motion between different sublattices. 

This, however, is not the only process that can take 
place. Fig. (|22|l shows that it is also possible for a hole 
to move two steps without creating a confining string 
potential. This is a second order process in perturba- 
tion theory. First we create a second hole two sites away 
from the first hole. We do this by moving that electron 
one site closer, next to the first hole. This gives rise to a 
double occupied intermediate state with a high energy U. 
Then, we let the same electron move again, now remov- 
ing the first hole. The amplitude for the total process 

.2 

is given by a = jj, much smaller than t, the first order 
hopping amplitude. However, the final state has exactly 
the same energy as the starting state. There are eight 
sites to which a single hole can hop in this fashion. Note 
that throughout the entire process, the magnetic quan- 
tum number in Fig. is preserved, e.g. S z ot = h/2 if 
perfect Neel order prevails everywhere around the frag- 
ment displayed in the figure for a lattice with an even 
number of sites. Within the t-J and t — J z approxima- 
tions to the Hubbard model, same sublattice hops are 
more masked. In the t-J model, only a third order pro- 
cess links a hole to one of its next nearest neighbors on 
the same sublattice. Albeit holding for models derived 
from the Hubbard model on bipartite lattices, in many 
instances the sublattice parity principle becomes much 
more alive and transparent within the original Hubbard 
model itself. 



Numerical simulations of the Hubbard and Heisenberg 
model with a single hole have found a discrete number 
of sharp peaks in the spectrum which could be identi- 
fied with string-states [25-43]. The longer the string, 
the higher the energy of the peak. Because there is a 
gap between the first and second peak, we can confine 
our attention to the lowest peak. This peak has a finite 
quasi-particle weight. The dispersion of this peak as a 
function of k, tells us about the effective movement of 
a hole through an antiferromagnet. A recent numerical 
result from Louis et. alS& for this dispersion relation in 
the Hubbard model is shown in Fig. I|23|) . 

The central result of numerous investigations is that 
they all found that a sharp quasiparticle peak appears at 
the bottom of a broad continuum of the hole spectrum. 
These quasiparticle poles form a coherent hole band with 
a width of order of 2J over a wide range of J/t, and 
the coherent propagation is made possible by "healing" 
a string of flipped spins by quantum fluctuations. 

One can easily Fourier transform the numerical low en- 
ergy hole dispersion relation e(k) to unveil the important 
real space quasiparticle motions and their respective am- 
plitudes. If we make the simplifying assumption that the 
hole is the quasiparticle then we find that essentially all 
of the low energy weight is distributed amongst the next 
nearest neighbor motions linking the hole to its sublat- 
tice. This is a consistent logical outcome of a strong local 
Neel order fused with the fact that the magnetic moment 
is a conserved quantum number. That the nearest states 
on the same sublattice (and those further away) have the 
highest weight could hardly be surprising. Hopping am- 
plitudes to sites which are further and further away on 
the same sublattice drops significantly with distance^. 

Stated alternatively, if on the two dimensional lattice, 
the hopping amplitude to any of the four colinear sites 
(twice (up, down, right or left)) is to and if motion to 
any of the four diagonal sites has amplitude £2) then the 
dispersion relation will read 

£fe = to — 4i2 cos k x cos k y — 2io(cos 2k x + cos 2k y ). (A2) 

From lowest order perturbation theorjii^, we imme- 
diately expect i2 = 2<o = — 0{\j) where the relative 
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FIG. 23: Quasiparticle band structure for a single hole on 
12 x 12 clusters of the square lattice with periodic boundary 
conditions and U = 8t. The solid line corresponds to the fitted 
dispersion relation (see text). The inset shows the bandwidth 
as a function of t 2 /U for U > 8t; the fitted straight line is 
— 0.022£ + 11. lit 2 /U. From Louis, Guinea, Lopez Sancho and 
Verges—. 



factor of two originates from the two paths by which we 
may reach diagonal sites by two consecutive hops as com- 
parde to the single two step route to longintudinal next 
nearest neighbors. As illustrated in Fig. I|24l) . renormal- 
ized hopping amplitudes of the same order of magnitude 
O(jj) reproduce the detailed and tedious numerical fits of 
Fig. (|23l) remarkably well. More generally, if one Fourier 
transforms the elaborate data encoded in Fig. (|23[) we 
find that the bulk of the Fourier weight corresponds to 
next to nearest neighbor motions with hopping ampli- 
tudes of the same order as anticipated from lowest order 
perturbation theorjiSS. 

In accordance with expectations from perturbation 
theory, it is indeed found numerically that diagonal 
(nodal) hopping is almost twice as large in amplitude 
as compared to longitudinal hopping (a property which 
might easily transcend to larger length scales by virtue 
of many of the small length scales in these materials) 
with both amplitudes of order 0(\j). In Fig. 
we show the theoretical dispersion curves coming from 
equation (|A2|I for t = -0.70|j and t 2 = -1.52|j and 
eq = — 68.15i. These curves capture all the essential 
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FIG. 24: Dispersion relation for a single hole in an anti- 
ferromagnet according to the simple view that the hole is 
a quasi-particle that moves on one sublattice by hopping in 
steps of two f equation IA2t . Here we employ eo = — 68.15i, 
ta = -1.52^ and to = -0.70^ with U = 8t. 



features of the numerical dispersion curve shown in Fig. 

It is obvious that the Brillouin zone has doubled as a 
result of the fact that the hole is moving in steps of two 
through the lattice. 

The very good fit obtained to the dispersion curve of 
a single hole in an antiferromagnet by introducing the 
three parameters (eo, to and £2) suggests that the com- 
plicated problem of a single hole in an antiferromagnet 
can be reduced to a single particle problem. This is a ma- 
jor simplification. If we really can neglect spin-flips and 
string-states, the single hole problem has been reduced to 
a single-particle problem. The spins have become static 
and their only function is to create a checkerboard like 
sublattice structure that keeps the hole moving on only 
one sublattice by forcing it to move in steps of two. Large 
J means in principle both large J z and J±. A large J± 
leads to spin-flips being relatively important and this will 
restore the string-states. Increasing the value of J, and 
thus of J±, leads to spin- flips which will restore a string 
to the normal antiferromagnet. However, the relevant 
region of the Hubbard model and the t-J model is the 
region of small values of J. 

Stated in more general terms, the constraint on the mo- 
tion of the holes (equivalent to the assumption of Neel 
order) may be regarded as a low energy selection rule on 
the matrix elements of the effective Hamiltonian. The 
total z component of spin (sublattice parity order) is a 
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conserved quantum number and transitions between dif- 
ferent sublattice states are banned. Whenever Neel order 
is strong, the states are effectively separated into two dis- 
joint Hilbert spaces each for a different sublattice parity 



(magnetic S z ) quantum number. As we demonstrate in 
the main text, stripe order may be interpreted as an im- 
mediate consequence of sublattice parity quantum num- 
bers. 
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